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Abstract . We investigate the effects of relatively rapid variations of the boundaries of an overmoded 
cavity on the stochastic properties of its interior acoustic or electromagnetic field. For quasi-static 
variations, this field can be represented as an ideal incoherent and statistically homogeneous isotropic 
random scalar or vector field, respectively. A physical model is constructed showing that the field 
dynamics can be characterized as a generalized diffusion process. The Langevin-Ito and Fokker— Planck 
equations are derived and their associated statistics and distributions for the complex analytic field, its 
magnitude and energy density are computed. The energy diffusion parameter is found to be proportional 
to the square of the ratio of the standard deviation of the source field to the characteristic time constant 
of the dynamic process, but is independent of the initial energy density, to first order. The energy drift 
vanishes in the asymptotic limit. The time-energy probability distribution is in general not separable, 
as a result of nonstationarity. A general solution of the Fokker-Planck equation is obtained in integral 
form, together with explicit closed-form solutions for several asymptotic cases. The findings extend 
known results on statistics and distributions of quasi-stationary ideal random fields (pure diffusions), 
which are retrieved as special cases. 
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1. Introduction 

Increasingly complex environments for propagation of acoustic or electromagnetic (EM) waves call 
for more powerful, efRcient and specialized techniques to characterize fields and associated quantities. 
Multiple scattering and diffraction, dynamic and complex geometries or configurations, multi-path 
propagation, wide-band quasi-random modulation, etc., all make deterministic analytical/computational 
or experimental characterization techniques increasingly less efficient. Moreover, in such applications, 
these conventional methods are often intrinsically inaccurate, because they typically presume single plane- 
wave excitation or localized idealized sources. This raises issues of analytical tractability or numerical 
stability, even in nonresonant "open" environments. As an alternative approach, a growing trend exists 
toward the use of statistical methods in such cases, e.g., for characterizing random scalar and vector fields 
inside overmoded resonant enclosures, as test or operational environments in their own right, or as field 
generators or simulators of complex and multi-scattering environments. A primary example is the mode- 
tuned or mode-stirred reverberation chamber (MT/MSRC) (e.g., [2]-[l]), which after an established 
career in acoustics is gaining increasing importance for application in electromagnetic compatibility 
(characterization of immunity, emissions, shielding, and absorption characteristics), high-intensity radio- 
frequency field generation, measurement of EM properties of materials, antenna characterization, or as a 
simulator of inhomogeneous time- varying propagation channels for radio waves in multi-terminal (MIMO) 
wireless communication systems [5]-[8]. In its basic configuration, a MT/MSRC is a highly overmoded 
resonant cavity whose boundary or excitation conditions are perturbed dynamically - thus generating a so- 
called "stirring" or mixing process - either mechanically (e.g., through continuous rotation of a reflective 
paddle wheel or other diffractor that is large compared to the wavelength and exhibits an acoustically or 
electromagnetically reflective surface; through rotation or vibration of cavity walls; etc.), electronically 
(by modulating the excitation signal in phase or frequency), via excitation by a noise source, or through 
any combination of the above. More generally, such a cavity can be employed for precise modelling and 
measurement of spatial and temporal distributions of EM energy [9]. An ideal MT/MSRC generates a 
statistically isotropic, homogeneous, incoherent, and (in the EM case) unpolarized field, which can be 
represented by an isotropic angular spectrum of random plane waves. In this way, it defines a canonical 
echoic EM environment that is the counterpart of unbounded (anechoic) free space. A characteristic 
feature of a MT/MSRC is the extreme sensitivity of its interior field to variations in the boundary 
conditions, at any location, akin to wave chaos displayed by non-integrable cavities ("billiards") with or 
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without time-varying closed or partially opened boundaries, e.g., |10|-[20j. As a result, mode stirring 
gives rise to a hybrid, i.e., amplitude-plus-frequency modulated interior field [21| . [22j . 

Fields in the presence of rotating scattering bodies or surfaces in unbounded environments have 
been intensively studied for several decades, in particular for rotations of a cylindrically symmetric 
and axially rotating scattering object, e.g., [13]-[2S]. More recent efforts have focused on irregular or 
complicated shapes inside enclosures. Pertinently, for an overmoded cavity in which a stirring process 
evolves arbitrarily slowly (adiabatically) relative to the effective intrinsic relaxation time of the cavity 
p6l[27] (that is, the weighted average of the modal decay times of all cavity modes that contribute when 
no stirring occurs), the overall system remains EM quasi-stationary. On the other hand, if the time 
scale of the stirring process is of the order of the intrinsic modal relaxation times (but not necessarily 
on a relativistic scale, and evolving still slowly compared to the period of wave oscillation), this may 
cause the interior field to "slip", i.e., to cause nonstationarity !22j. In this case, the source (excitation) 
field and the resulting interior (cavity) field exhibit rates of fluctuation that are of comparable order of 
magnitude, so that the imposed perturbations are no longer tracked faithfully (i.e., on the same time 
scale) by the induced interior field. In MT/MSRCs, this situation occurs most readily for acoustic waves, 
because of their significantly lower velocity of propagation compared to that of EM waves. As will be 
shown, the nature and parameters of the cavity field then become markedly different from those in the 
quasi-stationary case. 

In this paper, we study quasi-monochromatic random fields that are subjected to continuously time- 
varying boundaries of a highly overmoded resonant cavity, generating a nonstationary random field at any 
interior location. The framework is that of stochastic classical wave mechanics and stochastic differential 
equations (SDEs). Statistics of the field, its magnitude and energy density are obtained from their 
probability distribution functions (PDFs) as solutions of the associated Fokker-Planck equation (FPE). 
The energy density is of special importance because the incoherency of random fields makes this the 
fundamental EM quantity. The focus is on waves generated by a transmitting source (Tx) and perceived 
by a receiver (Rx) that are both located inside the cavity, as a canonical scenario for measurements by a 
field sensor, membrane, antenna, as well as radiation from, or susceptibility of an acoustic or electronic 
device under test. 

The paper is organized as follows. Based on a physical model developed in section [2l the SDE is 
derived in the present context in section[3l together with expressions for the drift and diffusion coefficients 
of the field in function of configurational and excitation parameters. These are employed to derive the 
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FPEs for the magnitude and energy density of the scalar and vector complex analytic fields in section |3] 
which, together with sectionlHl forms the core of the paper. The FPEs are solved in closed form, yielding 
the general PDF in integral form. The asymptotic PDFs for first-order (early) time dependence are given 
and their physical characteristics are discussed. The particular case of an asymptotically fully developed 
nonstationary process is analyzed in section [5] The results are extended in section [5] to second-order 
systems, including the case of first-order systems responding to a nonstationary field. The case of time- 
dependent but deterministic coefficients of the SDK and FPE is treated in section [T] Finally, a summary 
of the main results with conclusions is given in section [8l 

A selection of preliminary results for the case of a scalar field was summarized previously in [T]. 
The present paper follows on from a companion paper in which the nonlinearity and distortion of 
electromagnetic fields caused by a mode stirring process were studied from the point of view of random 
modulation of the interior field 22J. Experimental measured results demonstrating nonstationarity in a 
MTRC induced by accelerated motion of cavity boundaries have been reported in [55]. 

2. Transient mode-stirred fields 

The problem at hand is the temporal transition of the local interior cavity vector field Y_{t), from its 
initial stationary state at time t = to to its value at a next state at i = to + ^t- This transition takes 
place under the action of an external perturbation influencing an otherwise stationary cavity field X(t] 
and is governed by intrinsic characteristics (modal time constants) of the cavity and by the rate of change 
(velocity) of the action itself. The functional relation between Y{t) and X_{t) will be considered in section 

13 

If the field were responding instantaneously to changes in the boundary conditions - disregarding any 
pure propagation effect, i.e., phase delay - then the rate of fluctuation of Y{t) would be solely governed 
by that of the source (itself assumed to be stationary) at all times during the transition. In resonant 
environments, the situation is more complicated and rather intricate, particularly when coupling between 
modes exists, causing the transitions to become quasi-random. 

Assume XSt) and Y_[t) to be randomly amplitude- and hybrid-modulated time-harmonic fields, 
respectively |22| . Consider the transient field caused by an instantaneous step transition between two 
discrete boundary states taking place at t = to- At t = to-, the perceived field is in a well-defined state 
Y_{tQ-) whence it can be expressed as a weighted phasor sum of a finite number of vector components, 
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viz., the ampHtudes of the participating instantanous cavity eigenmodes [321 1301 or plane- wave vector 
components of a reahzation of the angular spectrum generated at Iq, respectively (cf., e.g., [STIIS^ V 

nto) = Y.ys*o). (1) 

2=1 

At any t, the local field can be represented as an analytic field with time-varying phasor Y_(t) in I-D or 
3-D complex vector space for scalar or vector fields, respectively, irrespective whether or not the field is 
quasi-stationary ^FH I30j . 

When, at t = to, a perturbation of the boundary conditions takes effect, each currently contributing 
mode or plane- wave component y^{t) decays, on average^, in accordance with its own characteristic time 
constant t^^'' governed by the cavity state at t — to. These prior contributions to Y_{to + At) will be 
further denoted with a superscript "(0)", as j/I^H^o + At). 

During the decay of each prior phasor, posterior phasors y^"'^' associated with the next-state 
(perturbed) cavity geometry start building up, each one governed by its next-state characteristic time 
constant T^^\to). In this process, energy is being exchanged between the participating modes owing to 
mode overlap in a lossy cavity [331 [M], in other words t^^-' ^ 0. This coupling may cause a phasor to 
have its amplitude temporarily increasing before on average decaying [35j . In other words, the decay of 
?/(°^ {t) or even F'"' (t) is typically not monotonic. The decay and build-up processes are complicated, 
because the eigenmodes effectively switch during these transients (birth-death- mutation processes), while 
their influence is prolonged as a result of finite rise and decay times owing to continuity of fields. Thus, 
at any posterior instance to + At, the transition is governed by sets of developing or fading modes: 

Y{to + At) = j2 + At) + E yf^ + At) 

i=i j=i 

= Y^''\tQ + At) + Y<^^\to + At). (2) 

We shall refer to this form of nonstationarity as stirring slip. It causes the phasors y^{t > to) to decrease 
exponentially, on average, each with different time constants while phasors undergo a relaxed and retarded 
change (they locally 'stretch', 'shrink', and rotate) as a consequence of variations of the cavity boundary. 

The associated phasor diagram is sketched in figure[T]for a scalar complex field Yq, = Y^—jY^. For a 
vector field, a denotes an arbitrary Cartesian component x, y or z, and l,^ represents a unit vector in the 

§ The details of the early decay of the effective field depend also on the exact instance (phase) of the y'*^' nt t = tg. Since 
this phase fluctuates randomly between recurrent events, each decay process exhibits largely different early decays for each 
participating phasor and, hence, for the resultant phasor sum between different starting times. 
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direction of that component. If the system remains quasi-stationary during the transition, ^^(io + ^t) 
is the sum Ya^^ (t) of a set of emerging posterior phasors (t) and Ya^^ (t) of a set of vanishing prior 
phasors yai{t) for to < t < to + At. For general nonstationary processes, however, Ya{to + At) is arrived 
at by a more comphcated dynamic superposition of a set of on-average decaying phasors {yai{t > to)} 
and growing new phasors {yajit > to)}- In case of "slowly" modulated signals, the situation can formally 
be described by an angular spectral plane-wave expansion for harmonizablc functions 

Mt) = / / cxp (-jfc • r) dX{t) (3) 
J Jn 

in which dX(t) is a process of nonorthogonal increments [36] . In the asymptotic limit of a pure Bachelier- 
Einstein- Wiener-Levy (BEWL) diffusion process (random walk), the y^a] [t > to) consist of independent 
increments building onto the existing prolonged {ya}{t > to)} = {Uaiito)}- 
Three particular regimes can be distinguished: 

• Quasi-stationary field: In this case, Ya{to + At) is statistically identical to ^^^(^o)- This can still be 
conceived as a special case of a random walk, viz., one that returns to the origin before each step: 
y^a-ito + At) = 0. Hence, Yi°\ta + At) = 0; cf. figure [TJi. In this case, 

Yaito + At)=Y^^\to + At). (4) 

• Fully developed nonstationary field: Here, Ya{t) is governed by ya}{to + At) — ya}{tQ) and can be 
represented as a traditional random walk within the complex plane, as shown in figure [1)3, whence 

Y^ito + At) - Y^'\to + At) + Y^"Hto). (5) 

In other words, Ya^\to + At) corresponds to the increment at t = to + ^t. 

• Intermediate nonstationary (slipping) field: This case is defined by y^a} {to + At) lying "in between" 
Uaiito) and 0, as depicted in figure [TJ:. Here, Ya^\to + At) "builds onto" a decayed version 
Y^°\to + At) of Y^°\to), whence 

Ya,{to + At) = yji)(to + At) + c„ri°)(to) (6) 

where Ca is a complex constant. 

Thus, only in the quasi-stationary case does Ya{to + At) not contain a contribution by Ya^\to + At). In 
any case, the discretized Ya{t) is a first-order Markov process because the transition between Ya{to) and 
Ya{to + At) is fully specified by Ya{to). 
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Figure 1. Models for transition between consecutive states of a random scalar wave in 
the complex plane: (a) quasi-stationary process, (b) fully developed nonstationary (BEWL) 
process, (c) general nonstationary process. The phasors y (to) represent individual principal 
components (modes, spectral plane- wave components, etc.) with random aplitudes and phases. 
In figure (c), the dotted phasors represent the corresponding phasors of the BEWL process of 
figure (b), shown for reference only. 

A description in terms of instantaneous cavity modes is only feasible insofar as the notion of 
instantaneous (eigen)frequency is applicable and sufficiently sharply localized. This, in itself, limits the 
applicability of the model to sufficiently small levels of nonstationarity [22 . These instantaneous modes 
are time-limited wavelets, rather than harmonic functions of infinite duration, enabling characterization 
in terms of an evolutionary spectrum [37j . 

To avoid the need for a parametric description 38, ch. 10], we shall assume that the field fluctuations 
are independent and much more rapid than those of the state of the system and its parameters. The 
latter can then be treated as constants. For the effective field, this is often a justifiable approximation. 
In section [71 we treat the case of time constants with arbitrary but deterministic time dependence. 
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3. Stochastic differential equations for real fields 

In this section, we consider the relationship of an output field Y{t) generated by an ideaHzed random 
input field X{t). Both X{t) and Y{t) are real modulated fields in the time domain; extension to complex 
analytic fields will be made in section ID 

3.1. Linear time-variant filtering of quasi- stationary random fields 

Consider a real scalar quasi-stationary fluctuating acoustic or EM field X{r,t) inside an enclosure, 
produced by a process of arbitrarily slow configurational changes (adiabatic boundary variations). This 
interior field X{t) is generated by a given time-harmonic EM source xo expQcjf) and satisfies the Hclmholz 
or vector wave equation with quasi-statically varying boundary conditions. The rate of fluctuation of X(t) 
is expressed by its correlation length Tp,x calculated from the autocorrelation function (ACE) px{t, At). 
The quasi-randomness of X{t) manifests itself only on time scales well in excess of Tp,x, by definition. 

The quasi-stationary field X(t) serves as input to a dynamic process of realistic configurational 
changes, i.e., occurring with nonzero velocities of the boundaries. The corresponding output is the 
resulting perceived "stirred" field Y{r^t) which is, to a certain extent, a weighted accumulation or 
aggregation of X{t). The transformation from X{r,t) to Y{r,t) is here restricted to be linear but 
possibly nonuniform, i.e., Y(r_, t) can be obtained from X(r, t) via linear time- variant filtering 39J. Here, 
we assume that this filtering is characterized, locally, instantaneously, and in the mean, by a first-order 
process with characteristic relaxation time constant r [40 . Assuming r 3> Tp.x, we can consider Y{t) to 
be a first-order Markov proces^U with an apparently white input process X{t) (cf. section l3T2|) . However, 
since X{t) is a physically realizable sample-continuous noise field, the Stratonovich picture is preferred 
for the stochastic formulation of the sample-continuous noise field. Eor completeness, however, we shall 
also list corresponding results for the Ito formalism [32] ■ 

In the limit Tp^xjT 0, the integrated noise B{r,t) ^ X{r,t')dt' represents a BEWL 
nonstationary process |43| . familiar from the theory of Brownian motion and random walks. This special 
case is analyzed in section [5l Eor practical purposes, X(t) and Y{t) may be considered as sample- 
continuous random processes, e.g., by considering an ideal white noise process subjected to a 'small' 
amount of local averaging [43l [44]. All integrated processes below are then properly defined in the 

II This is not a fundamental limitation, because in the other case a generalized Fokker-Planck equation can still be derived, 
cf. e.g. [m eqns. (1.27) and (4.1)]. 
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mean-square sense. 

We shall further be concerned with local statistical properties only, whence the dependence of X and 
y on r will be dropped, thereby assuming that all fields are evaluated pointwise without local averaging. 

3.2. Langevin-Ito equation 

Our first task is to derive the differential equation governing Y{t) with X{t) as its source term. To this 
end, we subdivide the observation interval [fg, tn] of length T = t„ — to into n equal subintervals [ti, ti+i] 
(i = 0, . . . , n — 1), each of length 5t = i^+i — t^, chosen such that 



The input field X{t) can be approximated by a series of discrete mutually independent Heaviside step 
functions and performs random jumps at discrete regular times ti. We assume that the levels X{ti) are 
maintained between ti and t^+i (sample-and-hold), i.e.. 



Since t ^ Y{t) \s & smoothed weighted aggregate of X{t). The fluctuations of Y{t) are a result of 
fluctuations of X{t) undergoing relaxed random spatial rotations and scalings of X(t). Consequently, 

Tp,Y > Tp,x. 



Tp,x ^5t<^T. 



(7) 



X{t, < t < t,+i) = X{t^). 



(8) 




-0.5 
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Figure 2. Discretized sample-and-hold input process X{t) and resulting output process y{t). The 



actual continuous- time stirring process is obtained as the limit St — ^ 0. 



At t — U^i^ two changes occur simultaneously (figure [2]): (i) the previous value of the input field 
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X{ti < t < ceases to exist. This can be interpreted as the "switch-off" of the source due to elapsing 
of the previous configuration. By itself, it would cause a decay in \Y{t > to the null field (ground 

state); (ii) the next value of the source field X{t > ti+i) comes into force. This can be interpreted 
as a "switch-on" of the source due to the generation of a new configuration. By itself, it would cause 
a transition of Y{t > i^+i) from zero toward the next equilibrium value X(ti-|_i) (asymptotic level if 
X{t > were to persist for t +00). Both effects are in competition if |X(ti-|_i)| > |^(i)|, while 

reinforcing each other if < |X(t)|. The resultant field is therefore 

'^^^ X{U). (9) 



Y{U<t<U+i)=Y{U)exp' * 



1 — exp 



In this model, the infiuence on the future state at ti+i stretches only as far back as the end value of its 

immediate predecessor state at ti, i.e., Y{t) is a first-order Markov process. Moreover, the field does not 

make discontinuous jumps between consecutive states and hence Y(t) is also sample-continuous. Defining 

dY{t) ^ Y{t) — Y{ti), ^ can be rewritten as the stochastic difference equation 
dY(t) exp(-^)-l l-cxp(-^) 

In the limit St ^ [implying Tp^x — * 0, on account of ([7])], (fTU]) yields the stochastic differential equation 

(sdeJs 

Y{t) = -T-^Yit) + T-^X{t). (11) 

This equation has been encoutered in a multitude of other applications, cf., e.g., [35 , but note the 
multiplicator 1/r for the source term here. Since r is independent of X{t) and, hence, of Y(t), the SDE 
pT|) is of the Langevin-Ito type and has solution 

Y{t)=yocxp(-*—^)+T-' f Xit')exp(-^—^)dt' (12) 



T 



where yo =^ Y {t ~ to) . The solution (|12p expresses the fact that Y{t) is an accumulation of all past values 
of X{t) weighted by an exponential relaxation. In certain scenarios, it may be possible to intervene 
manually and reset yo to zero. In general, however, Y{ti) equals the value of Y{t) that was reached at 
the end of [ti-i,ti] and is hence governed by Y{t — ti-i), St and r. 

3.3. Moments, drift and diffusion coefficients 

Apart from providing information on centrality, dispersion, etc., the moments of Y(t) enable 
determination of the coefficients of the Kramers-Moyal and Fokker-Planck equations for its probability 

*i If the process is irreversible 1151 . then it can be easily shown that ((TlJ generalizes to Y{t) = -(1/t(") )y (t) + (l/r'i) 
where r(0)'(i) = limt, _-t^_ ^^^^ rit). 
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density function (PDF) fY{y,t)- An ideal random, i.e., Gauss normally distributed source field X{t) can 
be completely characterized by its mean value {X{t)) and autocovariance function, with the aid of 

{X{t)) = (13) 
and, for an assumed wide-sense stationary X{t), 

{X{t)X{t'))^alpx{t-t'). (14) 
In the limit of ideal white noise, 

px{t-t')^5{t-t'). (15) 

In practice, the noise X{t) is a coloured [px{t) 7^ ^{t)]- In such cases, the idealization to white noise 
may be insufficiently accurate, whence the full expressions of px{t) must then be used. Several different 
approaches exist, in particular for colouration by an exponentially correlated X(t) |54]. Pure Markov 
processes associated with such an ACF are, however, not mean-square differentiable, as is physically 
required, because px (t) is not differentiable at r = and therefore require at least some small degree of 
local averaging [13], |44j . 

From (HH), for to <t,t' <ti, 



= yo exp 



(Yit)) 

(y(t)r(t')) = yo'exp 



r 

t + t' -2t 



0, 







Vo exp 



ft+t' ^2to 



t + t' - 2io 







hi: 









t — u + t' — v\ a\ 

—px[u - v)av 



exp 



*+oo) 



2t 



exp 



2t 

\t-t'\ 



\t-t'\ 



exp 



i + t' - 2io 



(16) 

du 

(17) 

(18) 

(19) 



where was obtained by considering the limit of a delta-correlated X(t). Hence, 



af (t) ^ {Y\i)) {Y{t)y = ^ 



1 — exp — 



2 {t - to) 



Recall that, in arriving at (|17p . the invoked interchange of the order of ensemble averaging and integration 
is permitted provided that a constant m exists such that \Y{t)\ < m and provided {Y{t)) exists for all t 
within the integration interval [5T1 section IV. 5]. For {t — to)/T ^ 1, the mean value and variance evolve 
approximately linearly with respect to t: 



{Yit)) ^ 2/0 1 - 



t-to 



<^Y — ^ {t - h) 



(21) 
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In practice, the stirred field Y{t) is typically a narrowband hybrid amplitude-plus-phase modulated 
function of time, whose statistics as a result of the stirring process vary slowly relative to the fluctuations 
of X{t) (weak nonstationarity) . Nevertheless, the application of the analytic field concept Y{t) does not 
restrict this process to be narrowband, for it suffices that the maximum frequency present in the spectral 
density does not exceed one half of the instantaneous average center frequency. 

The mean rate of variation of stir states (average derivative) of the stirred field is expressed by the 
drift coefficient -Dy''(j/,i), 



Y=y 



(22) 



in which y has a sharp, i.e., definite value at time t. The mean rate of increase of the uncertainty of the 
state (spread of derivative) defines the field diffusion coefficient Dy^y^t), 

[Yit + St)-{Yit))f) 



ot — >0 



which, unlike is independent of y. Since the source term X{t)/T in pT|) is independent of Y{t)^ the 

spurious drift {\ /2)dD^^\y ,t) / dy is zero. Consequently, the Ito and Stratonovich forms of the SDE (fTTj) 
are identical in this case |42j . 

Two particular asymptotic modes of operation can be retrieved as limiting cases of (ITT|) . If 
D^y^ = then we obtain a pure diffusion of the random field, because then Y{t) = t ^X{t) yields the 
BEWL process Y{t) = T^^B{t). At the other extreme, if r ^ (implying no local temporal averaging, 
i.e., zero memory) then Y{t) = 0, whence the original point process is retrieved, i.e., Y(t) = X{t). 



2! St 



Y=y 



'Y 

2r2 



(23) 



3.4. Fokker-Planck equation 



On limiting the Kramers-Moyal equation for finite differences, viz. 



5fY{v,t\V(i^h) 
St 



f-; n! \Sy 



{[Y{t + St)-Y{tr) 



St 



fYiy,t\yo,to) 



(24) 



to terms up to and including second order (n < 2), and upon substituting (|16p - (|17p . we obtain 



SfYiy,t\ya,ta) 
St 



l-exp(-f) 



St 



^[yfYiy,t\yo,to)] 



Sy 



2St 



1 — exp — 



St 



"_x_ 
2t 



1 — exp I 



25* ^ 



25t 



S'^fYiy,t\yo,to) 



Sy"^ 



(25) 
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In the limit (5t — > 0, this results in the FPE 

dfYiy,t\yo,ta) ^ _ _^ r 

dy 

- 

rdy 



dt 



yo,to) 



D'^\y,t)fY{y,t\yo,to)\ (26) 



2 2 



(27) 



dy"^ 

T 2t2 J dy 

with specified boundary conditions /y (±oo, i|2/0i ^o) = 0. This PDE governs the evolution of the transition 
PDF (TPDF) /y(±c3o,%o,io) or the PDF/y(y,i) whereby /y(y, t = to) = /y (2/0,^0), up to second order. 
The deterministic part of X{t) and hence of Y{t) satisfies the Liouville equation, i.e., (|27|) with cr^ = 
and r > 0, yielding a continuous probability current (cf. ([62])]. 

Using standard methods [ISl section 1.11], the general solution of (P7)l is obtained as an Ornstein- 
Uhlenbeck process 



fY{y,t\yo,to) 



1 T 


TTCT^ ^ exp 


2(t-to)' 
r 


} 



exp 





- 2/0 exp 


(-^)] 


2 




1 — exp 


2(t-to) 
r 


} 



(28) 



If, when considering a collective of boundary states, yo does not exhibit a sharp value across the 
ensemble but is characterized by a random variable Yq that itself exhibits the same stationary PDF, i.e., 
if 



■ exp 



X 



'X 



(29) 



then /y(2/,t) is itself stationary throughout, with 



exp 



D^Y^y^t) 



c 



+00 



fY{y,t\yo,to)fYo{yo,to)dyo 



■ exp 



'X 



Ty 
2 

X 



(30) 



as obtained from (j28p for {t — t^j/r +00, in which C is a normalization constant. This stationary 
solution ~ sometimes referred to as the subensemble extracted from /y(2/,t) as defined by j/o ~ exists 
provided the ergodicity condition for Y{t), i.e., 

l.i.m.T^+00 [hT{y) - {hriy))] = (31) 



lim 



+CXD r-\-oc 



— 00 'J —00 



WT jy) WT (2/0) 

2^2 



Jy {y, t\yo, to) - /y {y, t) /y (2/0, to)] di/dj/oj dMto = (32) 



is satisfied, in which wx[y{t)] = T^^ h[y{t)]dt is an arbitrary deterministic function. In particular, the 
difference between the moments for boundary averaging and spatial (interior domain) averaging [i.e., for 
WT^y) = 2/™] should fall off sufficiently rapidly with increasing sample length T. Thus, upon averaging 
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over all possible values of a Gauss normal Yq, the perceived variance is independent of t and inversely 
proportional to r: 

(af->=af-(^)=f . (33) 

Compared to fviy, t) for sharp initial values yo, the drift now vanishes in the mean and the field variance 
increases by an amount 2(7^ exp[— 2(t — to)/T], i.e., by a factor 2(7^/{exp [2{t — to)/T] — 1}. 

4. Application to analytic fields, energy density and magnitude 

Complex analytic scalar field 

Thus far, the stochastic characterization was for real harmonic or quasi-harmonic fields X{t) and Y{t). 
We now associate with X{t) the Gabor analytic field [22 ^{t) ~ iT^i^it)], again defined for real t, where 
Tl[X{t)] = n^^ X{t) / {t — u)du denotes the Hilbert transform of X{t). We shall further denote 

this complex field by X{t) ~ X'{t) —}X"{t), where X'{t) and X"{t) are the in-phase and quadrature 
components. Similar definitions and notations apply to the complex Y{t). 

If Xit) represents a vector field, its Cartesian components Xa{t) = X'^{t) — (o^ = 2/; ^) ^-^^ 

assumed to be mutually uncorrelated. Consider Xa[t) to be circular Gauss normal and delta-correlated. 
The output analytic field is then governed by two scalar SDEs, i.e., 

Y:}'Kt)+T-'Y:}'Hi)=T-^X'^'\t) (34) 

with {X'^'\t)) = 0, {Xa{t)Xa{t')) = a\J{t ~ t'). The solution is an Ornstcin-Uhlcnbcck process 

with the initial value ya'a Ya '\t = io). The associated mean value and covariance are 

(i^a^'HO) =2/;Wexpf-^V (36) 



f ( t + t'-2to\ , '^x^.c 
/ exp I 1 J — 

In the limit t/r, t'/r — > -l-oo 



f \t~t'\\ f t + t'-2t( 

exp — exp ' 



(37) 



2 

(y^W(t))--0, (y^(')(i)y^(')(i'))^!g:iexp(^-^). (38) 



4-. 2. Energy density 

4-2.1. Scalar fields or Cartesian components of a vector field 
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4-.2.1.1. Moments The energy density is the physical quantity of fundamental interest, because an ideal 
random field is incoherent, i.e., delta-correlated with respect to time. Consider the normalized energy 
density Ua{t) — Y^{t) + y"^'^(t), after normalization by a factor relating to the constitutive properties 
of the stationary medium, which are not important in the statistical characterization for a deterministic 
medium. The stationary Cartesian density decreases for increasing r: 



(2 ^lA 


2{t-to)' 






T 



(39) 



where = y'^^ + y^^^ = ^-nd a\ = a\i — cr\n. On account of the Isserlis moment theorem for 



a circular Gauss normal {t) , it follows that ( Ya 



3 K' 



, whence 



alJt)^2mt)f-4(Yi'{t))(Yf{t) 



(40) 



.2.1.2. SDE From 

{U^it + 5t) U^it)) = (Y^'{t + 5t) y^'(i)) + {Yfit + St) Yf{t)) 

2{t-to) 



Vo 



exp 



2t 

2 

y'ao ) exp 



2t 



T 

2(t-to) 



1 — exp 
1 — exp 



2St 

T 

2St 

T 



it follows that the energy drift coefficient is constant with respect to Uo 

i^«(..,t).limiM±i!)^^ 

" St^O ot 

where = Ua{to) = y^^. The energy diffusion coefficient is calculated in [Appendix A| as 



2 /^L A 


2(t-io)" 


= Wqo exp 






T 



Djj'(ua,t) = lim 



V^{t + 5t)-Uc.{t)X 



2St 



(41) 



(42) 



(43) 



Ua=Uc 



The SDE for Ua{t) can be obtained from and is nonlinear. Although the transformation from 

the complex analytic field to the real energy density is a multivariate one, the SDE can in this particular 
case be uniquely determined via inversion [HI section 3.4.1] because Y^{t) and Y^{t) are orthogonal, 
whence the diffusion dyadic D-Y^Y'i^ya^Va^^) diagonal. In the Ito formulation, the SDE is 



Uc.it) = D\jl{Ua,t) + ^ B^{t) 



'Xc 



^x, 
^ao ) exp 



2{t-to) 



(44) 
(45) 
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in which we have written Ba{t) for X^it), whereas in the Stratonovich formulation, 



idD 



(2) 



dUa 



exp 



4. 



2t 
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(46) 



(47) 



involving a multiplicative nonlinear white noise source process that is no longer of the Langevin-Ito type, 
which can be traced to the fact that D^j depends on Ua{t). Since the transformation from Ya{t) to Ua{t) 
is nonlinear, the Stratonovich formulation contains the additional spurious drift term — cr^ /t^. The 
SDEs (|45l) and ([47|) can easily be integrated numerically, for example using the Cauchy-Euler method, 
to generate sample functions. 



4.2.1.3. FPE Following the associated FPE for fu^{ua,t\uao,to) is 

2 /crl^ \ I 2{t-to) 



d 



2al 92 



with specified initial value Ua{to) = Uag, i.e., initial condition 



exp 



Uafu^ iUa,t\Uag,to)] 



(48) 



(49) 



/(7„ iUa,to) = S{Ua- Uao) 

and boundary conditions 

fu^{0,t\Uao,to) ^ fu^{+00,t\Uao,to) =0. (50) 

In ([35]), the product of the field variate with its PDF occurs through its second derivative, rather than 
its first derivative in (|27[) . The sought fu^ {ua , t) can be obtained via an integral transformation of 
(cf. [Appendix Bp and can be expressed as the Laplace inversion formula 

fua,{ua,t) = (j27r)"^ / Ts {s,t)exp(uas)ds 

J 7-joo 

for 7 > 0, where 

Ts (s,t) = exp 



(51) 



X exp 



'x^ 



'x^ 



exp 



exp 



{t-to~ t") 



"^xj 



dt" 
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X exp 



exp 



At'" \ dt' 



X exp 



1 A MqoT , 



(52) 



If the imposed boundary condition is (0+, t) = for all t, then J^s(s, i) is limited to the first term ii 
(j52p only. If, in addition, (wq, to) = for all Mq, then 



J^s{s,t) — exp exp 

E + oo [ 



- io) 



o-^ s 



Ei 



'2 , ^ T 

-{t- to) + 





T 




-Ei 




}) 









1 r 11^^^ 
f(t-to)-^jEi[-f(t-to)-^J} / 



exp 



El 



/n! 



(53) 



exp[-f(t-to)] E:ro 

where Ei (z < 0) = — J^^°°[exp(— z')/z']dz' is the exponential integral function 48, eqn. (8.211.1)] and 
El (z) = /+°°[exp(-z')/z']dz' = exp (-z) YI^o (-l)"n!/z"+i = -Ei (-z). With the aid of the residue 
theorem, ([5T1) -(|52 p can also be expressed as (cf. [Appendix B[ ) 

fu„ {ua,t) = ReSs=o [^s{s, t) exp (u^s)] 
corresponding to the coefficient of in the Laurent series expansion of Ts{s, t) exp (uqs), i.e., 



(54) 



+ 00 I 

J^s(s,t)exp(ucs) = UgS - exp 

1=0 



2 (t - to) 



E 



'x„ 



In 1 + ^ 7 



E 



■E 



n=l p=0 \ P 



2 (t - to) 



which contains infinite numbers of both positive and negative powers of . 



/£! (55) 



4.2.1.4- Limit PDFs In the hmit {t - to)/T 0, the exponential time dependence in (gS)), (jB.ip . (|B.5p 
and (jB.lOp disappears. The limit PDF can be obtained in closed form using a transformation of variables 
in the FPE (cf. [Appendix CP yielding 

{Ua + Uao) 



2 (t - to) 



Xc 



exp 



2 (t - to) fTx„ 



2 



(t-to)a2 



(56) 



valid for (t — to)/T ^ 0, which involves in general a Bessel function of negative fractional order 2i>, with 
v given by (jC.SP : 



< 0. 



(57) 



[Recall that /2iy(-) 7^ -^-2iy(-) unless 2;/ is an integer number.] Observe that is in general not separable 
with respect to Ua and t, unlike for quasi-stationary fields. 
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As a special case, if u^a'T l<y\ ^ 1 in addition, then /i = 1/2, = 0, whence ([55]) reduces to a PDF 



for a diffusing and drifting energy density characterized by a Nakagami-Rice distribution with PDF 



2 (t - to) (T^^ 



■ exp 



{t - to) a 



(58) 



2 (i - io) t'x^ 

whose variance increases hnearly with time. The order (type) of the PDF has asymptotically increased 
to zero. This PDF constitutes a special case of the general solution ([?T|) . For arbitrary t, it corresponds 
to a noncentral PDF including a constant bias Uao- Note that for Ua„ — 0, (|58|) reduces to a negative 
exponential (^2) PDF for which /c/„(0+,i) 7^ 0, i.e., it does not strictly satisfy the prescribed boundary 
condition at Uao — 0. Therefore, the latter PDF represents a singular case. 

At the other extreme, if (t — to)/T —>■ +00 (limit of long time or vanishingly short memory) then (|48p 
becomes independent of Uq^, viz., 

This FPE coincides with (jC.2p if we choose /i = 0,j^ = — l/2in (|C.5|) . An asymptotically stationary 
solution exists, because the coefficients of the FPE (|48| become asymptotically independent of time. 
Hence, from (1551) . the PDF for {t ~ to)/T ^ +00 approaches 



fu^ (Uq,<) = 



■ exp 



2{t- to) <tI 



h 



(^ - to) 



(60) 



2 (t - to) 

Comparing (|C.2[) with ([55]) . and (|56p with it is found that the change in drift gives rise to a 

different functional dependence of fu^ (m^, t) on u^. Thus, as time progresses, not only do the parameters 
and moments of the distribution evolve, in a continuous manner, but also the nature of the PDF itself. 
Both features are discussed in more detail below. In the special case when v = —1/2, i.e., Ua^ = cr^^/r, 
the initial and final asymptotic PDFs coincide. 

At the other extreme, the asymptotically stationary state [t/r +00) corresponds to Ya{t) 
performing a random walk in the complex plane with full return to its reference state (uaoi^o) after 
each completed step (complete instantaneous relaxation). The underlying PDF of the diffusing Ya{t) is 
then of course identical to the one for Xa{t), viz., a circular Gauss normal PDF with time- independent 
parameters. Therefore, (155)) reduces to a central x| PDF, as expected: 



■ exp 



'x^ 



In general, the probability current for fu^{ua,t), i.e. 



]uAUa-,t) = DljHua,t)fu^iUa,t) - 



U, 

4. 



dUa 



1 d 

2 dua 
2u„ 



D^^l{Ua,t) fuS'^a,t) 



'-fu^iUa,t) 



(61) 

(62) 
(63) 
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is the sum of a convection current iiafuaWcnt) plus a diffusion current dfu^{ua,t)/dua of probability 
as the PDF evolves in time. Expressions for ju^ follow with the aid of the above expressions for fu^ , but 
are in general rather cumbersome. For the special case Uaa — 0, however, 



■ exp 



2{t - to) V 2(i - to)a|^ 
which vanishes in the limit {t — to) / t —t +00. 



1 



2(f-to)a£ 



(64) 



4-2.1.5. First-order time- dependence of PDF: asymptotic PDF Having investigated the solution 
/c/q {ua,t) of ((48|) for asymptotically small and large times, we now consider the dynamics of the PDF as 
a function of time. As mentioned before, the PDF for general t is obtained by inversion of (|52p . or ([55)1 
in particular. With the boundary condition /;7„(0+,<) — 0, the kernel ([52)1 reads 



Ts (s,t) = exp 



X exp ■ 



2{t-to)crl 




{'-'if:) 



t" 



■ exp 



--[t-to 



2cr^ s 



-t' 



dt" 



(65) 



/itl|+]] 



For < (i - to)/T < 1, on defining t'" ^ - r2/(2cr^^s) and witltLI t"'/T < 1, the exponential function 
in the kernel can be expanded to first order in {t ~to)/T with reference to t'" — 0, yielding 



Ts{s,t) 





2{t~ta)' 




j exp 


L 




T 






-dt" 



2(t-to) " 



in which t' = 2aji {t - to)/T^ [cf. (ICTl) ] and 



-1 



2{t-to) 2{t~to)uo 



(66) 
(67) 

(68) 

+ Note that this must be valid for all values of s when integrating the inverse Laplace transformation l|51|l . In particular, 
for the part of the Bromwich contour parallel to the imaginary axis in the complex plane, this requires an appropriate limit 
process for the principal value of the integral when — > 0. The contribution of the essential isolated singular point s = 
in the inversion integral can be neglected if it is assumed that /t and the positive SRfs] = 7 in USU I are both sufficiently 
large whence (a\ /t)7 > 1 and that 2t"' /t < 1 for t'" G [0,t - to]- 
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Figure 3. Evolution of short-time asymptotic probability density functions of the energy density 
for scalar field {p = 1) . 



Except for an additional s-independent scaling factor and the fact that the exponent of (1 + si') is now 
time-dependent, ([ST]) is of the same form as (|C.3p . Consequently, as an extension of we now obtain 



exp 



valid for {t — to) / t <ti 1, in which the time-dependent parameter 



2{t- to) a 



[t " to) a 



(69) 



v{t) = 



(70) 



exhibits a linear dependence on (t — to)/T. Again, j/(i) becomes constant with respect to time when 



(71) 

7 

This corresponds to f in ([57]) taking on the value —1/2 and, hence, (|69p like ([5^ to coincide with the 
limit PDF ([50)1 at all (early) times, as indicated before. 

Figure [3] shows the asymptotic PDF for selected values of {t — to)/T <C 1 with <o = 0, r = 0.3 
s, crx^ — 0.3 V/m. The corresponding dependencies of the mean, standard deviation, and coefficient of 
variation are shown in figure 3) 

When T/a\^ 1, the term t/((7^^s) in the exponent of (1 + st') in ((66l) must be accounted for. 
On application of contour integration in the complex s-plane [cf. [Appendix E| with the definitions (|E.3p - 
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Figure 4. Evolution of instantaneous mean value, standard deviation and coefficient of variation 
of the energy density for scalar field {p = 1). 



E.5I1 ]. fu^{ua,t) can be expressed as (|E.24p . i.e., 



2(t-to)a 



2{t - to)a 



Xo, 



X sin < vr 1 — 



x-l 
^ 2{t-to) 



exp {—Uax) dx. (72) 



4-2.1.6. TPDF for random initial value (extraction of subensemble) For a xi'distributed random initial 
energy density Uag exhibiting the same statistics as Ua itself, i.e., 



fu^o (wao) = exp — 

^Xa V ^Xa 

the PDF /[/^(uQ,,i) coincides with the asymptotic xl-distributed TPDF: 

p + CO 

fu^iUa,t) = / fu^{Ua,t\Uao,to)fu^giUao)'iUao 



(73) 



ai{l + exp[-i^]-exp[-^]} 



X exp 



TUo 



ai{l + exp [-^^]-exp [-'-^]} 
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(t-to)<CT 



exp 



TUo 



(74) 



which, unhke (|30p . remains time-dependent. 

For more general random the issue of estimation of its PDF arises. This can be achieved in 
a systematic manner using a Bayesian scheme, i.e., starting from an assumed prior PDF /(/^^(uao) ^^'^ 
adding information of the observed /(7q (wq, t|uQ,o, io) to yield an improved posterior estimate: 

lu^aiUaaWa) = C fu^{Ua\Uao) fu^g {Uao) (75) 

where C is a normalization constant. This expression can be used in an iterative scheme to update the 
posterior distribution. For the prior distribution, one can use a known or a chosen ad hoc (Gauss normal 
or other) distribution centered around a classical (e.g., maximmn likelihood) estimate of the mean value. 



4.2.2. Vector field In the following, we shall use the subscript t to explicitly denote quantities associated 
with the total (i.e., vector) field, consistent with the use of the subscript a before. For consistency, results 
will be expressed in terms of statistics relating to Xt rather than Xa, in order to yield self-sufficient 
distributions and statistics. 

For the energy density of the mixed total field lt(i), we readily find in analogy with sections 14.2.1. II 
and l4.2. 1.21 that 

mt)) = (Y^'it)) + (Y:"it)) - ^ ^ (76) 

n(i)^„, . 
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(77) 



(78) 



where cr^ = cr^ 



g\ = 3(Tv and Ut = Ux + Uu + Uz = 3Ua. The associated SDE follows as 



'X, - -"^x, 

4. 



f^t(i) = -< ( ^-«t„ )exp 



in the Ito formulation, or 
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2 {t - to) 
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Ut{t) 



6t 



Utit) 



Bt{t) 



(79) 



(80) 



in the Stratonovich formulation. The corresponding FPE is 
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dt 
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'3^0^ [^t/t/t {ut,t\uto,to)] 



(81) 
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with uto = Utito), i.e., 

fut{ut,ta) ^S{ut- uto) 

and 

fuAO,tWto,to) = fUt{+00,t\uta,to) = 0. 
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(82) 



(83) 



Its transformation to a first-order ODE and its solution via Laplace transformation are obtained in a 
similarly way as in [Appendix B[ but some subtle changes in the coefficients arise. Details are given in 
[Appendix D| and lead to the general solution again being given by the inversion formula (j51[) but now 
with 
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(84) 



If /(7t(0, t) = then Ts{s, t) is limited to the first term in (|84p only. If, in addition, fut{ut, to) = then 



J^s{s,t) = exp < —3 exp 



(t - to) 



3r 



Ei 



- (t - to) + — 



Ei 
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(85) 



In either case, /;7t(ut,t) is obtained via inverse transformation of (|84p . or (|85p in particular. 

As a special case, we consider again the early-time or large-relaxation limit {t — to)/T ^ 0. Upon 
scaling t as 

t'--^(t-to) 



3t^ 



the FPE (ED becomes 
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whence ([M)) simplifies to 

Ts {s,t') = exp 
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1 + st' 

UtaS 



exp 



-3 1 



3i» 



expl-^j(l + .t') . 

Except for the exponent of the second factor, ([55]) is identical to the result for the scalar case (jC.3[) 

Therefore, we can again use (jC.4p with a = t' ^ + s,/3 = y/u^/t' , v = /i — (1/2), but now defining 

3Mt„r 



1/ = 1 



With this choice of parameters, the limit PDF of Ut for {t — to)/T^Q is 

3t2 f utY r 3t2 (ut + uto)l r r3T^yMrwt;; 
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Although this result is formally identical with ([55]) . recall the different definition of v in 
to (jC.5|) . If, in addition, UtgT/a'j^^ <C 1 then /i = 3/2, i^ = l whence (|M| specializes to 
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(90) 
compared 

(91) 



which is to be compared with (|58p . By contrast, the limit {t — to)/T +oo yields the same limit PDF 

as the one for Ua, i.e., (j60p . after replacing with wtp, viz., 
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4-2.2.1. First-order time- dependence of PDF For {t — to)/T ^ 1, the evolution of the PDF for the total 
energy density can be analyzed in a manner analogous to that for the Cartesian energy density in section 
14.2.1.51 As mentioned before, the PDF for general t is obtained by inversion of ([52|) . or ([53|) in particular. 
Instead of ([65l) - ([67|) . we now have with fu^ (0+, t) = that 

UtoS 
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in which t' is again defined by 

1 

The resuhing PDF shows again a time-dependent order, i.e., ([90]) generahzes to 



and 



where 
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3ut„T 3{t-to)uto t-to 
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If r/cr^ ^ 1, then ^ generahzes to [of. (|E?24l) with the definitions (|EJ)) - (|R9)) ]. 
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Figure [5] shows the asymptotic PDF (f97|) for selected values oi {t — tQ)/T <^ 1 with to = 0, t = 0.3 s, 
(Tx„ = 0.3 V/m, ~ 0.9(7^^. The corresponding dependencies of the mean, standard deviation, and 
coefficient of variation are shown in figure [6l 



4-2.3. Generalized field On occasion, it is expedient to handle energy densities with more than three 
degrees of freedom, for example when collating data for i statistically independent locations or for 
differently generated realizations. There are then in general 2p = 2m£ degrees of freedom, where m = 1 
or 3 for single-point Cartesian or vector fields, respectively. The analysis and results of section [4.2.2l can 
be repeated mutatus mutandis. Denoting Up{t — to) = Up„, the Stratonovich SDE and the FPE are now 
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and 
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(101) 



pT^ dup 

respectively, with fu-p{up,to) = 5{up — Up^^) and with the boundary values fup{0,t\upg,to) = 
/[/p(+cxD, t|up„, to) — 0, and the p-dimensional circular Gauss normal white noise process Bp defined 

by 



Bp{t) 



\Bpm 



ELi mt)]' + mt)Y 



(102) 
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Figure 5. Evolution of short-time asymptotic probability density functions of the energy density 
for vector field (p = 3) . 



in which Bn = B'^ — j-B". Corresponding to ([69|) . we have the solution 



now with 



pr 



2(t-to)(T2 



exp 
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(103) 



(104) 



4.2.4- Extraction of subensemble An ideal random initial energy density Uo = \Yq\'^ exhibits a X2~ 
distribution with standard deviation 2(7^. Performing an additional averaging, the Up^ in the above 
expressions is to be replaced by (Upg) = 2cr^ /''' = P'^x /'''■ This coincides with (|7T|) . whence the type 
of fij^ (although not its statistics) then becomes again time-independent. Hence we retrieve (j60|) with 
replaced as indicated. For the total or generalized field, however, the order of the Bessel function 
in fij^ remains time-dependent and the PDF does not approach (j92p when [t — t^j/r ^ +oo, except if 
<j\ /{pt) = (^"x^/t in which case v{t) — (p/2) — 1. 
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Figure 6. Evolution of instantaneous mean value, standard deviation and coefficient of variation 
of the energy density for vector field (p = 3) . 

J^.S. Field magnitude (envelope) 

The instantaneous field magnitude (amplitude) A{t) of an analytic process Y{t) is defined by the envelope 



y'^(t) and can be obtained from U{t) via the variate transformation = ^JU{t). 



4-3.1. Scalar fields or Cartesian components From (|42p and (|A.6p . we obtain the drift and diffusion 
coefRcients of Aa{t) as 



D^ia^.t) 
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fdA^ 
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(105) 
(106) 
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in which Oa = \ya\ \/u^. The SDE for Aa{t) is thus quasi- linear. In the Ito formalism, this equation 
is obtained on application of the Ito differentiation formula [49] : 

dU^{t) (dC/„)2 



dA40 = ^dt + ^dC/„ + i?-^ (dC/„)^ 



at 
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in which pS]) is substituted for and with {dUa)"^ — {iUaCrj^ /r^)di, yielding 



'-Ya 2 

a^. 



exp 



2(t~io) 



2r 



(109) 



This SDE can also be obtained directly as Aa (t) = D]^ 



2D^^^ /ax^ ) Ba{t). A form of (fTOQ]) solely 



involving amplitude (rather than field) statistics can be obtained by replacing ax^ with cr^^/-\/2 — (ti'/S) 
for a Rayleigh distributed [43]. In the Stratonovich formulation, the last term in (|105p is omitted. 
Substitution of (gT]) into Aa{t) = f/Q(i)/[2Aa(i)] then yields again Thus, unlike for [/^(i), no 

spurious drift of Aa{t) occurs, as is also apparent from the fact that Dj^ is independent of a^. 



For the special case of a process dominated by relaxation (t +oo), (|109p reduces to a result that 
has been obtained previously using a different approach [52l eqn. (4.4.41b)]. 
Using pUS)) and (dUTD, the FPE for Aa{t) is obtained as 
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dt 
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/Ac {o-a,t\ O-ao i ^o) 

2r^ dal'^- ^^^^^ 
As an alternative to solving (|110p . the (T)PDF for /a^ can be determined more straightforwardly from 
variate transformation of the corresponding distribution for Ua- For example, the limit distribution of 
Aa{t) when {t — to)/T and a^j^r/cr^ <C 1 is obtained from ([55]) as 



/a„ (oq, i|aQo, io) 
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2 (t - to) 
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, (oq, > 0,t > to). 



(Ill) 



This distribution has also been encountered in a related problem of one-dimensional power dissipation in 
electronic circuits subjected to external illumination by random fields [50j . 



4.3.2. Vector field For At{t) we obtain, using ([77]), ([75]) and the definitions (|105p and (|107p applied to 
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The corresponding SDE is 
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whose short-term fundamental solution, assuming an unbiased underlying field, is given by 
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4-3.3. Generalized field For a generalized field with 2p degrees of freedom, the SDE and FPE follow by 
straightforward extension as 
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(117) 
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2pt2 ' 

respectively, with short-term fundamental solution, again assuming an unbiased underlying field, given 

by 



/Ap {ap,t) = 



pap 



(t - to) a\ 

5. Example: BEWL diffusion process 
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exp 
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(119) 



The above results for the dissipative Ornstein-Uhlenbeck diffusion can be specialized to the case of a 
Bachelier-Einstein- Wiener-Levy (BEWL) process Bp{t), i.e., a pure diffusion, corresponding to the limit 
t/r in the above (cf. figure [1)3). In this case takes, the field requires an arbitrarily long time to relax 
to its asymptotic equilibrium state. Here, we formulate the results immedately in terms of the generalized 
field, from which local scalar {p = 1) and vector {p = 3) fields follow as particular cases. 



5.1. Energy density 

5.1.1. SDE For a BEWL process, we have that 



{Bp{t)) = 0, {Bp{t ~ to)Bp{t' ~ to)) ^^{t + f -2to-\t- t'\ 
Following an analysis similar to that in section 14.2.11 we obtain 

{Up{t + 5t) - Up{t)) = '^rJ^-upA St, 



(120) 



(121) 
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{[Up{t + 5t) - Up{t)Y) = — Up{t) St + Ml+p){St)^ + O [iStf] . (122) 

Combined with and applied to a generalized field, the SDEs for Up in the Ito and Stratonovich 
formulations become, respectively. 



Mt) ^K^-Upo + J ^BAt) ) , (123) 
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1 \ , IUp{t) 



(124) 



5.J.^. FPE Following (gB]), the associated FPE for Jup{up,t\upg,tQ) is 



J,p 



2 



K/c/p ("p,*l"po:*0)] (125) 

now with the boundary conditions fupi^i^WpoT^o) = /!7p(+oOi^l''^po7^o) = 0- 

With the former condition, it is verified that the probability current (j62p . which on account of current 
continuity can also be expressed as 
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- [u'ju, (u;)] du;, (126) 



p J 

vanishes at itp = for all t. This may be expected on physical grounds, on account of the positivity of 
energy [fupiup < 0,t) ~ 0]. In general, however, jUp{up > 0, t) 7^ except for the trivial case Up„ — 0. 

5.2. Field magnitude 

Through the variate transformation for the field magnitude Ap{t) cx yJUp{t), we obtain the drift and 
diffusion coefficients as 



^a:(«P'^) = ^ (128) 
whence the SDE in the Ito formulation is 



^,(,) . D% + ^B„m . (1 - -L) 4^ + ^ (129) 

Alternatively, p29p may be obtained from the Ito formula. The SDE in the Stratonovich formulation 
(2) 

is identical, because D^^ is independent of Up. The corresponding FPE follows again by substitution of 
(fml) and (fT28)) . 
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In certain practical applications, the goal is to determine how a (for sake of simplicity) first-order 
system with characteristic time constant t responds to a nonstationary cavity field. The analysis below 
complements a similar analysis for nonlinearity and distortion given in section IV. B of [22j . To this end, 
in section we first determine the response of a second-order system to an ideal white noise source field 
(quasi-stationary interior field), by converting a single second-order SDE into a system of two coupled 
first-order SDEs with well-separated time scales. This approach is akin to a widely used technique for 
analyzing a coloured noise process [54] . The result is then used in section 16.21 to find the response of a 
first-order system to a nonstationary effective interior field. Throughout this section, we assume a scalar 
field for simplicity, thereby omitting its subscript a for brevity. 



6.1. Quasi- stationary cavity field and second-order receiver system 

The results of sections |3] and |4] for a simple first-order Rx can be extended to higher-order systems. For 
a second-order system represented by an RLC equivalent circuit, the response is now governed by 

Y{t) + 2CY{t) -I- uolY{t) = v{t) + T-^B{t) (130) 

where v{t) is a deterministic excitation (e.g., direct illumination or biasing field), lo^ = (LC)^^ and 
2C — uio/Qo, in which Qo — lOqRC or [lUoRC)"'^ for a parallel or series RLC circuit model, respectively. 
The mean, variance, and covariance oiY{t) and Y{t) can be derived by converting the single second-order 
SDE (|130p into a system of two coupled first-order SDEs, as shown in [Appendix F[ yielding 
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(133) 



+ (C - v/C^^) {l - exp [-2 (C - VC^^) } 



o, ,2 
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(134) 
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87-4 (^2 



{l-exp [-2 (C + v/C^^) t] } 



+ {l - exp [-2 (C - \/e - •^i) } - 2 [1 - exp (-2^)]] . (135) 

Note that the mean value of the transient (ahhough not the variance) now depends not only on j/o but 
also on ?)o, i.e., the initial velocity of the mixed field. Thus, the result will in general depend on whether 
the reset occurs at the beginning of the stirring process (reset from "stand-still" ) , or during the process 
("flying-start" reset). As a generalization, for a third- or higher-order Rx a dependence on the initial 
field acceleration j/o or higher-order time-derivatives occurs, respectively. 

From (|R2|l . the FPE for the joint PDF /y ^(y, y) with v{t) = follows as 
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The stationary joint PDF fyy-iy, y, +oo) is obtained from (|136p as 



/y,F(y,2/,+oo) = /y(y,+oo)/i'(2/, +oo) = C'cxp —y exp - 
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i.e., jyy factorizes into two Gauss normal marginal PDFs for Y and Y . The corresponding FPE for the 
marginal PDF fY{y,t) is well known to be (cf. eq. (1.4.245) in [38 ) 
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The stationary marginal PDF /y(?/, +oo) can also be directly obtained from (|138p and is similar to the 
corresponding PDF for a first-order process, viz.. 



/y(2/, +Oo) = 



^ — exp 
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(139) 
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whence — > +00) = 0, (Tyii +00) = Qd {^lo'^T^) ax- As may be expected, lowering the Q- 
factor of the Rx is seen to result in a reduced variability of the perceived field. This effect, which must be 
accounted for, is irrespective of any loading effect that insertion of the Rx may impose on the stationary 
cavity field which excites it. The marginal PDF +00) follows from (jl37p and (|f 39p as 



whence /iy(t +00) = 0, cry(t — > +00) = Qo/ {2ijJoT^) ax ■ The latter distribution is relevant to the 
characterization of the rate of fluctuation, whose sample maximum value is important in applications 
[22j . The distributions of the energy density and field magnitude are qualitatively the same as for the 
first-order Rx. 

6.2. Nonstationary cavity field with first- order Rx 

The previous analysis presumed X{t) to be quasi-stationary, i.e., to represent an input field that is 
sufficiently slowly varying relative to the intrinsic time constant r^h of the mechanically static propagation 
environment (cavity). If the mode-stirring process is sufficiently rapid (e.g., for electronic stirring), 
however, then the field presented as input to a detector or test device with its own characteristic time 
constant r may itself become nonstationary. In effect, the nonstationary cavity field Y{t) now becomes 
the new input field X(t) for this device. This complicates an analytical treatment. A time-domain 
numerical calculation then seems the most appropriate approach. Alternatively, a complex frequency 
a +iLu may be introduced to incorporate transients or nonstationarity into this X{t) |22| . In order to gain 
further insight, however, we shall here pursue an analytical solution under the assumption that X{t) is 
short-term stationary (i.e., quasi-stationary within a sufficiently narrow time interval). Thus, in a certain 
sense, the cavity now "generates" its own transfer function for the stirring process as a result of the finite 
response time Teh of the chamber (cavity) which is long in real time compared to the scale of fluctuation 
Tp.x of the stirring process. 
Specifically, in this case. 




(140) 




(141) 



thereby assuming that 



Tp,x < Tch < r < 1. 



(142) 
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Each one of the first two inequalities in (|142p enables either one of the equations in p4ip still to be 
treated as a SDE driven by a white noise process (i.e., a source function that is fluctuating very rapidly 
relative to the response function). In (|14ip . X{t) — B{t) is the actual and possibly nonstationary cavity 
field, as would be perceived by an idealized instantaneously responding dot sensor or component [r = 0, 
i.e., Y{t) = X{t)\, while B{t) is defined as the time-integral of a now fictitious quasi- stationary cavity 
field Z{t). Equation (|I4ip holds in the mean, when the higher-order fluctuations of the field decay are 
neglected, and can be written as the single second-order SDE 

Y{t) + ^-±^Y{t) + (TTch)"' Y{t) = (TTeh)"' Bit). (143) 
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From (fn51) and ([^5]) with 2C 



+rT/ and uj^ = 



(rTeh)-\ the joint 2-D PDF /^^(y,?/) 



and marginal 1-D PDF fY{y,t) with the initial condition /y(yoiio) = <^{y ^ Uo) have associated EPFs 
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dt T dy ^ " ' 2t^ 

where the latter equation is valid to first order in the ratio (tcIi/t). The stationary joint PDF follows 

from (|f44p as 
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The stationary marginal PDFs follow from (|145p - (jl46p as 
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where we made use of the fact that Tch t. 

The general nonstationary solution can be obtained as a series expansion, but its expression is more 
cumbersome. We limit ourselves here to calculation of the time dependence of the first- and second-order 
moments. Following from (jF.ll|l and (jF.12p . we obtain 
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Since Tch <C r, the Rx experiences only a second-order asymptotic effect, viz., 



+oo) ~ (j/o + Tchyo) exp ( -- 
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with (fT55ll and pSTjl leading to the PDF (fT47| . The variances ([TSTl) and (fT58| indicate that nonstationary 
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mode stirring yields a lower root-mean-square (rms) magnitude of fluctuation (Y'^{t)) but a larger rms 



rate of fluctuation y (Y'^it)) compared to (quasi-)stationary mixing [i.e., stirring or tuning {tcIi/t 0)]. 
In the limit Tch — +oo, we retrieve the characteristics of a Markov process, i.e., ay +oo, as may be 
expected. Furthermore, (|155p and p56p show that no asymptotic bias occurs (/iy, /iy — > 0), although 
the initial bias (for i/r — > 0) now increases when Tchyo 7^ 0. 



7. SDE and FPE with time-dependent deterministic coefficients 



In the above, we assumed that the decay constant r was independent of the stir state at time t. For 
example, in the case of a nonstationary response of an electronic circuit immersed in a stationary mode- 
tuned field X{t), T is entirely determined by the circuit itself and is then indeed independent of t. By 
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Time 1 (units 51) 



Figure 7. Time dependence of the amplitude of the discretized nonstationary real-valued 
stirring process with state-dependent decay constant T(t). Dotted, solid and dashed curves 
represent the contribution to Y{t) when having reached steady-state at the starting time, i.e., 
Y{to — iSt) — X{iSt), i — 1,2,3. These contributions evolve with different time constants 
governed by the subinterval [ti,U+i]. 



contrast, when analyzing a mode-stirred cavity field resulting from nonstationary mechanical or electronic 
stirring, the chamber's decay constant t = Q/uj oc u:^^!'^ is governed by the particular set of participating 
eigenmodes at each stir state. These modes and, hence, r then depend on t (cf. figure [7]). In a first 
approximation for this scenario, we may assume r (being an average value taken across a typically large 
number of modes within a frequency band that is narrow relative to the center frequency) to be varying 
relatively slowly compared to the local N(t) and hence for r to be a quasi-deterministic but time- varying 
function of t (Bachelier process). Repeating the foregoing analysis of a discretized stirring process, 
under the often plausible assumption that the decay of |y(to-)| toward zero is governed by the same 
characteristic time as the competing growth from zero to |y(to- + At)|, and subsequently taking the 
limit (5< ^ 0, the SDE becomes 



Y{i:) = -T-\t)Y{t)+r-\t)X(t) 



now with the solution 
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which can be integrated for specified deterministic T{t). For example, assume that the rate of change of 
T is independent of t and denoted by = f, and that T{t) can be hnearized when taking the hmit 5t — > 0, 
so that we can write 



T{t) = {t-h)ui2+Tl 



for ti <t<t2, whence 
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For given T{t), the drift and diffusion coefficients can be explicitly calculated. In general, the resulting 
FPE is still of the form 



dfup{up,t) ,,.dfu^{up,t) d^fu^{up,t) 



(165) 



which can again be solved by Laplace transformation and conversion to a first-order ODE, by proceeding 
as in [Appendix B[ Instead of (jB.Sp , the appropriate transformations of variables are now 

s' = s + w{s,t), t'=t + u{s,t) (166) 

where we define 

wis, t) = g{to)sh - s2 /■ g(t')dt' (167) 

J to 

u{s,t) = [g{t)s]-\ (168) 

with which the coefficient oi dJ-s'{s' ,t')/ds' in (jB.lOp is constant with respect to t' and the coefficient of 
dJ-s'{s',t)/dt' vanishes again. 

In the case when the time dependence of r(t) is weak, we can immediately state the FPE in its 
first-order approximation, i.e., (|165p with 
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Nonstationarity of scalar (acoustic) and vector (electromagnetic) random classical fields in complex 
propagation environments has been modelled as a diffusion process and analyzed statistically based on 
a description in terms of stochastic differential equations, ft was shown that closed- form PDF solutions 
of the FPE for the energy density can be derived if the underlying field is ideal Gauss normal, without 
recourse to a typical polynomial expansion of the solution for the associated Sturm-Liouville problem. 
It was found that, in the general case, PDFs of the time-varying energy density are not separable with 
respect to the energy density and time variables. Limiting and asymptotic forms of these PDFs in the 
short- and long-time limits have been obtained. The main results are contained in ([TT]) . (|10ip ~ (|103p . 
pT7)) - pl^ . pM 1) - P^ . P^ - plD)) . The energy distribution of a generalized field was found to exhibit 
a Bessel / distribution with time- varying order and statistics. The results are, in particular, important for 
the characterization of sample statistics for nonstationary fields, whose random fluctuations are governed 
by dynamics that exhibit significant departures from those for stationary systems |55j . 
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Diffusion coefficient for scalar or Cartesian enery density 



In this Appendix, we derive the expression ([^5]) for the diffusion coefficient D^^\ua,t) for Ua{t) = 
^Q^(0 +^q'^(^) when Ya is circular Gauss normal at any t. To this end, we expand {[Ua{t + 5t) — C/q(<)]^) 
with respect to Ya'\t)- Making use of the fact that the process has independent increments 
Ya{t + 5t) — Ya{t) that are conditional to Ya{t) and, hence, to Ua{t), we obtain after straightfowarded 
calculation 

{[u^it + St) - u^{t)f) = (Kit + 6t) Y^{t)t) + (Kit + St) y^'(t)]^) 

+ 2 (K(t + St) Y^it)f) (Kit + St) Y^it)f) 

+ 4 [r^^ (0 (Kit + St) Y;,it)f) + Yfit) (K{t + st) - Y^'it)f)'_ 
+ 8 wr^W (Kit + St) Yiit)] Kit + st) - Y^it)]) . (a.i) 

The first pair of terms in (jA.ip involves higher-order moments that are in general difficult to express 
in terms of Ya '\t), unless the Yait) are circular Gauss normal. In this case, using the Isserlis moment 
theorem, we have that 



r^W(t + <5t)-y^W(t) 

For the third and fourth terms in (|A.1[) . 



(A.2) 



Y;('\t + st)-Y;('Ht)] ) =([y^W(t + 5i)] ) + {[Y;yHt)] y 2 (Y;}'\t + st)Y^^'Ht)) 



St 

exp 1 -- 



2r 



fcto 



exp 



2(i-to) 



-1 + 2exp 



St 

T 



exp 



25t 

T 



(A.3) 



In the limit St — ^ 0, to first order in St^ the expression (|A.3[) reduces to iu\ /T'^)St which is independent 
of y'ao to this order. Finally, the last term in (jA.ip vanishes, because Y^it) and Y^it) are uncorrelated. 
Thus, provided Y^^it) and Y^it) are identically distributed, (jA.ip becomes 

[u^it + St) - Uo.it)f) = 8 (Kit + St) r^W]')' + 4 + St) - u^it) (a.4) 

Uo,it)St + 0[iSt)^] (A.5) 



whence the diffusion coefficient for the Cartesian energy density is 



Djj'iua,t) = lim 



Vo.it + St) - Ua.it)Y 



2St 



(A.6) 



This coefficient is independent of Ua^ because the dependence of (|A.3|) on St is only of second order. 
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Solution of FPE for Cartesian energy density 



In this Appendix, we derive the general solution ([?T|) - ([5^ of the FPE-based boundary value problem 
O-dSOl), i.e., 



d 



exp 



2{t- to) 



du'i 



(B.l) 



for the Cartesian energy density Ua with initial value Ua{to) = Uag, i.e., initial condition 

fUa. {Ua,to) = 5 {Ua - Uao) (B.2) 

and boundary conditions 

fUa,{0,t\Uaa,'to) = /c/„(+00,^|Wao;*o) = 0. (B.3) 

Since Ua takes positive values only, we perform a Laplace transformation of the stochastic PDE with 
respect to Ua, viz., 

J^s{s,t) = / fu^{ua,t)exp{-Uas)dua, 3?[s]>7o. (B.4) 

JQ 

Here, 70 denotes the abscissa of convergence for fu^{ua,t), with 70 < +00 because fir^(ua,t) is of 
exponential order owing to the upper boundary condition in (|B.3I) . The transformation of (jB.ip is 



ds 



+ 



r2 dTsis^t) , 1 



2a 



dt 



'x^ 



exp 



2(t~ to) 



now subject to the corresponding initial condition 
J^S (s, <o+) = exp i-UagS) . 



'x^ 



exp 



T 

2{t- to) 



Ts{s,t) 

fuAo+,t) 



(B.5) 



(B.6) 



Specifically, by the physical nature of energy density, we are interested in solutions that satisfy the 
Dirichlet boundary conditions 

fu^ (0+, t) = 0, fu^ (+(X3, t)^0 (B.7) 

in accordance with similar properties of known PDFs for the Cartesian energy density of imperfect quasi- 
stationary random fieldj*] [30]) [IZ]. The second condition in (|B.7|) is the natural boundary condition 

* The former condition is not satisfied for the negative exponential PDF for the Cartesian energy density of a perfect 
reverberation field, for which /(7o(0+,t) = 1. For an imperfect field, more elaborate models |30l . 1471 show that the 
boundary condition is fulfillod, which underlines the generality of the solution found in this section in all practical cases. 
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for all PDFs; the former condition makes the right hand side of (|B.5p vanish (homogeneous PDE). The 

PDF will also satisfy the Neumann condition dfi/^{ua — Q+,t) /dua — 0. However, (jB.Sp is generally 
independent of its specific value dfu^{ua = 0+,t)/i9uQ, as a result of the pre-factor Ua for the term in 
d^fujduo:' in jRl]). 

The first-order PDE (|B.5P may be converted to an ODE of same order, through a suitable 
transformation of coordinates. Defining 



t' = t-U 



so that 



ds 



ds 



2(j\y dt' ' 



dt 



dt' 



(B., 



(B.9) 



the terms containing dJ-'s{s,t')/dt' in (|B.5|) then cancel, whence this PDE reduces to 



ds' 



-1 

with the re-scaled initial condition 

^2 



1 - 



T aj^ s' I 



2t' 

exp I - — 



exp(— Mqqs') = exp 



'x^ 



'2aj^J' 



.fuAO+,t) 



(B.IO) 



(B.ll) 



in accordance with (IB.6P and (jB.Sp . This yields the general solution 



J^5'(s',i')exp 



exp 



2ct^ t' 



-1 



2t' 

exp 1 - — 



ds 



X exp 



1 - 



1 

„/// 



'Xo, 



f 2t' 
exp 



'x^ 



/c/JO+,t) 



'x^ 



exp 



2t/_ 

T 



'x^ 



ds" 



ds". 



(B.12) 



A further transformation of (|B.12[) via i"^') = l{2a\ s"*-'^), followed by a transformation back to the 
original s and i, yields the final result given by ([5^ . 

The integral expression of the solution can be re-expressed with the aid of residue calculus. The 
point s = is an essential isolated singularity in the complex s-plane, because s~^ appears in the 
argument of exponential functions in J-'s{s,t). Therefore, calculation of (jSip using the residue theorem 
requires calculation of the Laurent series for J^5(s, t) exp (mq,s). The asymptotic points 7±joo are also 
essential isolated singularities, because this series contains terms with positive as well as negative powers. 
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on account of the factor exp(uc,s). The Bromwich contour has its vertical straight segment located at 
3fi[s] = 7 > and closure is across the negative (left) half of the s-plane. When tracking this contour, the 
transition across the positive real axis is governed by lim|s_^|_,o exp(l/s) = exp(l/7), which is defined 
provided 7 7^ 0; otherwise lim|s|_^o 6xp(l/s) = lim„-^+oo exp(— jn) which is indeterminate. The integral 
across the closing arc of the contour can be verified to satisfy the Jordan lemma for sufficiently small 7, 
whence its contribution vanishes. The inverse transformation then yields (|54p - (j55p . 
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Limit PDFs 



In this Appendix, we derive the early-time limit PDF (|56p of the Cartesian energy density governed by 
the FPE (gSl), i.e., ((RT|1 . For (i-to)/r 0, the exponential time dependence in (|RT|) . (|R5)) and (|RT0)1 
disappears. Upon scaling the time variable to 



t' ^ 



it - to) , 



the FPE (jB.ip transforms as 
dt' 

and ([5^ becomes 

Ts (s,t') = exp 

= exp 



S 

1 + St' 

1 + st' 



dfuAno.,t') ^ d^fuAuc.,t') 

+ Ua- 



dUa 



exp 



- 1- 



dul 



(C.l) 



(C.2) 



(1 + St') 



We now make use of the identity [ISl eq. (6.643.2)] 



T{^JL + v 



+ 00 



^ exp {—ax) (2f3y/x) dx 



(C.3) 



(C.4) 



r (21/* + 1) y z,u; / i~> \ '-'^ / JQ 

where Af_^i_i/(-) and l2i^{-) are Whittaker and modified Bessel functions of the first kind and order 2z/, 
respectively, in which we define the parameters as a = t' +s, /3 = y/Uao/t', and 

1 ^ Ua„T 

2 ^ 2o-|^ 



< 0. 



(C.5) 



With this choice, we have that M-^^v{z) = z^exp(z/2) with which the left member of (jC.4|l can be 



rewritten as {uao)'^t' exp{uao/[t'{l + st')]}/{l + st 



Hence (jC.4p becomes 



exp 



l + st' 



(1 + sf) 



exp(-» 

(Uag)'^ t' Jq 



(ua)" exp - (t' ^ + Ua hu ^ ^ ^" "" ^ dUa (C.6) 

which corresponds to (IC.3|) and yields the inverse transform of !Fs{s,t'). After a final transformation 
back to the original time variable, the PDF of Ua for {t — to)/T ^ is obtained as ((56)) . 
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PDF for energy density of vector field 



In this Appendix, we derive the kernel ([5^ for the general solution of the FPE-based boundary value 
problem ([5T|) - ([55)) for the energy density Ut{t) of a vector field Yt{t). Compared to the scalar case (jB.sp . 
the transformed FPE is now 

2{t-to)' 



ds 2a\ s2 dt s \ a 



) exp 

Xt 



T 

2{t- to) 



with transformed initial condition and Dirichlet boundary conditions 

J^S (s,io+) = exp (-utos) , 



(D.l) 



(D.2) 



/c/,(0+,t) = 0, /c/,(+oo,t) = 0. 
The pertinent transformation of variables is 



s' = s 
t' = t-t 







2ai- s 

with which (jP.ip - (jP.SP transform to 



dTs'is'^t') , 3 



ds' 



- 1 



exp 



-1+3 1- 



2£ 

T 



2t' 

exp I 



/[/.(0+,t) 



3t 



with 



3t" 
2a\^s 



exp(-utos') = exp 



2a\t' 



(D.3) 



(D.4) 



Similar to the procedure in [Appendix B[ the transformation i''^'^ = "ir"^ / {2a\^s"^'^), followed 
transformation of Ts' (s', t') back to a form in function of the original s and t, yields 



(D.5) 

(D.6) 

by a 
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Mellin integral for asymptotic PDF of energy density 



In this Appendix, we seek to calculate the inverse Laplace transform of 



F(s) = cxp(--^)(l + a.) 

1 + as 



via the Mellin integral 



fu (u) = (j27r)"^ / F{s) exp (us) ds, 7 = di[s] > 70 
J-f-ioo 



(E.l) 



(E.2) 
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where 70 is the abscissa of convergence, in connection with the early-time asymptotic PDF (I72p for the 
scalar or Cartesian field, and (|99p for the vector field. For the case of scalar or Cartesian energy density, 

(E.3) 



2(i - <o)ffL 



1 



1 



c = 1- 



'Xc 



(E.4) 

(E.5) 
(E.6) 



Mo = "ao 

with < {t ~ to) / T <^ 1, Uao > and < UaoT/a\^ <C 1, whence a>0, 0<5~1, 0<c<C I/mq;,, and 



< c/b <C I/mqq. For the energy density of a vector field, 

2(i-toK, 



a 



3t2 



b =i\l- 
c =3(1- 



3^itoT^ 
-I. ) 


[1- 








k 



2{t-to) 



Mo = Utg. 



(E.7) 
(E.8) 

(E.9) 
(E.IO) 



with corresponding conditions holding. 

With regard to the general case (|E.ip . the issue centers around the integration along the incisions Ci^± 
which surround and exclude the singularity — 1/a (see below). On account of the Casorati-Weierstrass 
theorem, exp[— mos/(1 + as)] takes on any given nonzero complex value in the neighbourhood of this 
singularity upon integrating. As part of a product of exponentials in F{s) exp{uas), however, this factor 
does not cancel when considering mirrorred locations on pairs of corresponding segments or arcs as part 
of the overall contour of integration. 

The special case for which mo = 0, as well as being of considerable practical interest, is easier to 
solve, i.e.. 



F{s) = (l + as) 



; exp 



In (1 + as) 



a, 6, c > 



(E.ll) 



which we shall consider first. Since — c/s] may take values smaller than one, the order condition 

F{s) = 0(8-''), \s\ +00, k > 1 (E.12) 

that would guarantee the existence of (jE.2[) . is not satisfied. Since this condition is merely sufficient for 
the existence of (|E.2p . an a posteriori check C[fi/^{ua)]{s) — F{s) for fu^iua) will therefore be required. 
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Li 


L2 


L3 


u 


6*1 = arg(s - si) 


IT 








— TT 


6*2 arg(s - S2) 




TT 


— TT 


— TT 


B3 = arg(s - S3) 


n 


TT 


— TT 


— TT 


60=^-62 + 63 


TT 


TT 


TT 


TT 



Table El. Argument of s £ C with respect to critical points Si when moving along the branch 
cut. 



To invert (|E.11[) . we apply contour integration in the complex s-plane (figure [ETJ. To this end, 
we extend F{s) analyticalljlj across the left half of the plane. Since {b — c/s) is in general (complex) 
fractional, we apply a branch cut to maintain a single-valued function. We choose this cut along the 
negative real axis, with the phase reference 6 = arg(s) = along the positive real axis: 

s e Cc = C \ {s| .s = X < 0}, arg(s = X > 0) = 0, - tt < arg(s) < tt. (E.13) 

Define the following special points from (jE.lip : 



S2 = 0, S3 



(E.14) 



where si is the critical point of (|E.11[) . S2 is the branch point, and S3 is a regular auxiliary point introduced 
to assist in the further analysis. Rewriting (jE.lip as 

F{s) = [|l + as|cxp(j0i)]l'-*l^'^P(j''°) 



= |l + asr-*l'="P^^''''^exp 



exp (j6'o) 



(E.15) 
(E.16) 



where, on account of a, 6 > 0, 

6*1 = arg(l + as) = arg(s-si) (E.17) 

60 =TT + arg - ^) = TT - arg(s - S2) + arg(s - S3), (E.18) 

Table IeT] shows the phase progressions when moving along the contour C around the branch cut, in the 
positive direction indicated in the figure. 

The joint contribution by the segments Li and L4 is 

|1 + ax\ exp (j7r)]l''~~l '''^P^^'^-' exp (uax) dx 



F{s) exp (wqs) ds 



lim 

Pi— >0,-R^oo 



-(l/a)-pi 



For notational simplicity, we shall denote this extension by the same symbol F{s). 
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x=Re[s] 



Figure El. Contour of integration C in complex s-plane. 



(l/a)-pi 



[|l + aa;|exp(-j7r)]l 



oxpQtt) 



lim { j2 

pi^O,i?,->oo 



-(l/a)-pi 



11 + ax\ 



exp (uax) dx j 



= j(-l)i-^+t2 r^^l + axyi"-^ 

A similar calculation for the other pair of segments i2, can be shown, with the aid of table IeI| to 
yield a pairwise vanishing contribution to the Mellin integral. 

Upon application of the Jordan lemma, the contribution across the arc Cr vanishes in the limit 
R +00, on account of 



sm 



TT {h ^ exp(MQx)dx. (E.19) 



lim (supc. |F(s)|) = lim (supp [1 + ai?exp (j6')] ^+ R<=^v(.m \ 



lim 
0, 



[l + ai^expQe)]^ 



exp 



i?exp(j0) 



(E.20) 
(E.21) 



because | exp(j0)| ^ for all 9. 

For the contribution along the circle C2 with radius p2, on substituting ,s = S2 + P2 exp(j0), we have 

lim [(s - S2)F(s)exp(uQs)] = ^ lim / F[s(p2)] exp [uas(/92)] dp2 = 0, (E.22) 
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on account of the integral limit theorem for vanishingly small p2 across an arc of given opening angle |53| . 
Note that S2 is a regular point [F{s2) — exp(ac)]. Likewise for the semi-circles Ci^±, using de I'Hopital's 
rule on the right member of (jE.ll[) . we obtain 



lim [(s — si) i^(s) exp (mqs)] =a 

s — ys-i 



-b-c 



' exp 



ft- to 



X lim exp 



cos 9 \ T 
In (pi) ft-to 



a 



exp (j6l) 
exp (j^) + Uapi exp (j6') 



(E.23) 



cos 9 \ T 

Although we have been unable to calculate the latter limit analytically, numerical computation for selected 
values oi 9 ^ ±tt/2 indicates that this limit is zero because of pairwise cancelling contributions. Hence 
the contributions by the arcs Ci,± vanish as well. 

Collating the above results and applying the Cauchy integral theorem on the closed contour C, we 

thus obtain a real integral expression for fuaWa)'- from (IE.2P and (|E.19p . 

(-l)i-^H 



{1 — ax) ^ sin tt ( H — J exp (— Uqx) dx 
i/a L V xJ\ 



(E.24) 



with which ^ and ([99]) follow by substituting (|K3)) - ((E31) or (|K7)) - ((R9|) . respectively 



Appendix F. Mean, variance, and covariance of quasi-stationary cavity field and 

for a second-order system 

In this Appendix, we derive expressions (|13ip - (|135p for the first- and second-order statistics for a general 
second-order system with a quasi-stationary input field, as characterized by (|130p . incorporating the 
particular case of a nonstationary cavity field filtered by a first-order system as specified by (|14ip . If we 
define 



Y{t) 
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-1 




' i = 






Y{t) 




[^l 2CJ 




















v{t) 




r-2 





(F.l) 



then (|130p can be rewritten as the Ist-order vector ODE 

tit) -g • lit) + hit) + gB{t). 
To integrate (jF.2p . we multiply by exp (ai) = J2n=o i—^)" /"■■' ^^^-^ 

exp (gt) ■ dVit) + exp {gt) ■ g ■ Y{t)dt = exp (gt) ■ [h{t)dt + gdB{t)] 



(F.2) 



(F.3) 
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Since the left hand side of (|F.3p represents d [exp (ai) • Y{t)\ , integration of this equation yields 

exp {gt) ■ Y_{t) = Y„+ [ exp (gt') ■ h (t') dt' + [ exp {gt') ■ gdB {t') . (F.4) 

J to J to 

The second integral in (|F.4p can be integrated by parts as 

I exp {gt') ■ gdB {t') = exp (gt) ■ g [B{t) - B{to)] - f B {t') d[exp {gt') ■ g\ 

J to J to 

= exp {gt) ■ g [B{t) - B{to)] - g- exp {gt') ■ gB {t') dt' (F.5) 

J to 

whence substitution into (IF.4|) and taking i3(io) =0 yields 

Y_{t) = exp {-gt) ■ Yo + gB{t) + [ exp [~g {t - t')] ■ [h {t') + gB {t')] dt' . (F.6) 

J to 

The expansion of the matrix exponentials (into an infinite series of matrix powers) can be avoided by 
making use of the matrix Laplace transform C [exp (— ai)] (s) = (s/ + a) ^ i.e., 

f +00 



exp {-at) exp (-st) dt = (s^ + 2(s + u;^ 



S + 2C 1 

: ,2 



(F.7) 



Taking the scalar inverse transform, element by element, yields 

exp (-gt) = (Ai - A2)"^ 

— Ai exp (Ai<) + A2 exp (A2t) + [exp {\\t) — exp (A2i)] exp (Aii) — exp {\2t) 

-Lul [exp (Alt) - exp (A2t)] -Ai exp (Aii) + A2 exp {\2t) 



(F.; 



where A1.2 = — C T ■\/C^ ^ ^-''^ the roots of + 2C,s + loI — 0. 

Alternatively, one may orthogonalize Y_. To this end, we seek an orthogonal transformation y' = c-Y , 
i.e., such that 

Ai 

A2 

Here, a primed quantity signifies a transformed (new) variable. This transformation is specified by 



c • a • c ^ = 



(F.9) 





Y'{t) 




A2 -1 




Y(t) 


(F.IO) 




- (A2-A1)-' 






Y'{t) 




_ -Ai 1 _ 




Y{t) 




in which Ai_2 as defined above are the solutions of det(a 


- A/) ^ 0. 


The associated vector FPE for Y^{t) 



yields a bivariate Gauss normal PDF with a mean value 



yo exp (Alt) 


, where 






yo 


= c ■ 




_ ?)oexp(A2t) _ 











(F.ll) 
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and covariance matrix 



where 



exp + \j)t] — 1 



'X 



4T4 [ul - C2) 



l-cxp(2Ait) 
2Ai 



l-cxp[(Ai + A2)f] 
A1+A2 



l-cxp[(Ai+A2)t] 
A1 + A2 



l~cxp(2A2t) 
2A2 



(F.12) 



(F.13) 



In terms of the individual processes Y{t) and y{t)^ the PDF /y(j/, t|y^, io) is also bivariate Gauss normal, 
with = c^^ • and ct^ = cr^ ■ a^, ■ {c^^)^ ■ This leads to expressions (|13ip - (|135p for the mean, 
variance and covariance of Y{t) and Y{t). 
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